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ON THE KERNEL FUNCTION OF THE INTEGRAL EQUATION RELATING LIFT AND DOWNWASH 
DISTRIBUTIONS OF OSCILLATING WINGS IN SUPERSONIC FLOW 1 

By Charles E. .'Watkins and Julian H. Berman 


SUMMARY 

This report treats the kernel junction of the integral equation 
that relates a known or prescribed downwash distribution to an 
unknown lift distribution for harmonically oscillating wings in 
supersonic How. The treatment is essentially an extension to 
supersonic How of the treatment given in NAOA Report 1284 
jor subsonic How. For the supersonic case the kernel junction 
is derived by use oj a suitable jorm oj acoustic doublet potential 
which employs a cutoff or Heaviside unit junction. The kernel 
junctions are reduced to jorms that can be accurately evaluated 
by considering the junctions in two parts: a part in which the 
singularities are isolated and analytically expressed, and a 
nonsingular part which can be tabulated. 

The kernel is treated jor the two-dimensional case, and it is 
shown that the two-dimensional kernel leads to known lift 
distributions jor both steady and oscillating two-dimensional 
wings. The kernel junction jor three-dimensional supersonic 
flow is reduced to the sonic case and is shown to agree with 
results obtained jor the sonic case in NAOA Report 1884, a/nd 
the downwash junctions associated with “horseshoe" vortices in 
supersonic flow are discussed and expressions are derived. 

INTRODUCTION 

In reference 1 the kernel function of an integral equation 
relating a known or prescribed downwash distribution to an 
unknown lift distribution for a harmonically oscillating finite 
wing of arbitrary plan form was treated for compressible 
subsonic flow. The purpose of the present report is to 
extend this treatment of the kernel function to supersonic 
flow. 

The kernel functions under consideration arise when 
lineorized-boundary-value problems for obtaining aerody- 
namic forces on oscillating wings are reduced to integral 
equations involving the distribution of pressure or wing 
loading as the unknown. In such integral equations the 
kernel functions play the important role of aerodynamic 
influence functions in that they give the normal induced 
velocity or downwash at any one point in the plane of the 
wing duo to a unit pressure loading at any other point in the 
plane of the wing. 

As the kernel functions arise in the analysis, they are 
mathematically defined by rather intricate improper inte- 
grals and possess singularities as high as second order. It is 
therefore desirable to isolate the singularities and determine 


their explicit nature in order to make the integral equation 
more amenable to solution, in particular amenable to solution 
by approximate or numerical procedures. 

Approximate lifting-surface theories for finite wings, such 
as the methods developed by Falkner and Multhopp (refs. 
2 and 3) and others, have afforded considerable success in 
the calculation of aerodynamic coefficients for steady sub- 
sonic aerodynamics. Similar approximate methods have been 
successfully employed to obtain coefficients for two-dimen- 
sional oscillating wings in subsonic (compressible) flow (for 
example, refs. 4 and 5) and are now being extended to the 
finite oscillating wing in subsonic flow by Harry L. Runyan 
and Donald S. Woolston of the Langley Aeronautical Labora- 
tory and by W. P. Jones (ref. 6). It is reasonable to expect 
that these methods can be further extended to apply to 
finite wings in supersonic flow. 

In supersonic flow, solutions of the boundary-value prob- 
lem for some particular plan forms and downwash conditions 
can be obtained in the form of infinite series in terms of a 
parameter involving the frequency of oscillation (see, for 
example, refs. 7 to 10) or in the form of rather complicated 
definite integrals (refs. 11 and 12). The infinite-series 
method furnishes a relatively simple means of obtaining the 
loading on oscillating ■wings for low values of the frequency 
parameter, but for large values of this parameter the series 
expansions converge so slowly that recourse must be had 
to other procedures for obtaining the wing loading. One 
feasible method is to study and develop approximate pro- 
cedures for solving the integral equations that involve the 
unknown loading and its associated kernel function. The 
first step toward such a development is to isolate and deter- 
mine the explicit nature of the singularities of the kernel 
function; this step is accomplished in the present report. 

The report contains the derivation of the kernel function 
in the form of an improper integral and a reduction of this 
integral to proper form. The singularities of the kernel 
function are isolated and expressed analytically, and the 
nonsingular parts are reduced to a form readily amenable 
to numerical evaluation, as was done in reference 1 for sub- 
sonic flow. Some expanded forms of the kernel function 
are derived, and one of these is used to obtain a reduction to 
two-dimensional flow. In appendix A, the limiting case for 
sonic flow is derived and shown to agree with the results in 
reference 1. Appendix B is devoted to certain integrals of 


i Supersedes NAOA Technical Note 3438 by Charles E. Watfclns and Julian H. Berman, 1958. 
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the kernel function. These integrals relate to “horseshoe” 
vortices in supersonic flow, as treated, for example, in the 
steady case by Schlichting in reference 13 , and may be of 
interest in certain modes of application. 


SYMBOLS 


c 

velocity of sound 

Jo, Ji(r) 

modified Bessel functions of first kind 

J»(x) 

Bessel function of first kind 

ICo, Ki 

modified Bessel functions of second kind 

K(x 0 , y 0 ) 

kernel function for three-dimensional flow 

Jc 

reduced-frequency parameter, fa/V 

K(xo) 

kernel function for two-dimensional flow 

v), m 

lift distributions 

La, L\ 

modified Struve function of first order 

l 

unit length 

M 

Mach number, V/c 

P 

r=Vi/o*+2 1 

perturbation pressure 

S 

region of xy-plane occupied by wing 

t 

time 

U(x) 

unit function 

V 

forward velocity of wing 

w(x,y,t) 

downwash velocity, e iu ‘w(x,y) 

w(x,y) 

complex amplitude function of prescribed ver- 
tical velocity 

x,V,z 

Xo=X— £ 

y 0 =y—v 

1 

Cartesian coordinates attached to wing mov- 
ing in negative x-direction 

5 (x) 

Dirac delta function 

Z,v 

Cartesian coordinates used to represent space 
location of doublets 

p 

fluid density 

<t(x,y,z,t) 

velocity potential, e ial ~4>{x,y,z) 

!>{x,y,z) 

complex amplitude function of velocity po- 
tential 

i(x,y,z,t) 

acceleration potential, e l “‘f(x,y,z) 

Kx,y,z) 

complex amplitude function of acceleration 
potential 

ei e 
II 

circular frequency of oscillation 


ANALYSIS 

INTEGRAL EQUATION RELATING DOWNWASH AND LIFT DISTRIBUTION 

The linearized-boundary-value problem for the determina- 
tion of the aerodynamic forces on a wing can be immediately 
reduced to a problem of solving an integral equation that 
relates downwash and lift distribution. The purpose of 
this section is to introduce and briefly discuss this equation. 

Since the integral equation has the same formal appear- 
ance for subsonic and supersonic flow and is derived in various 
publications (for example, refs. 1 and 14 ), the equation will 
not be rederived here but will be formally stated so as to 


serve as a starting point in the analysis. In keeping with 
linear theory, the wing is considered as a plane, impenetrable 
surface S which lies nearly in the xy-plane as indicated in 
sketch (a). 



Sketch (a) 


The x,y,z coordinate system and the surface S are assumed 
to move in the negative x-direction at a uniform velocity V. 

In terms of these coordinates, the integral equation may 
be formally written as 

w(x,y)=~- p JJ £(£, v)K(x-Z,y— v )dt dr, (1) 
a 

where w(x,y) is the complex amplitude function of prescribed 
vertical velocity or downwash at points in S and is defined as 
follows: 

w(x,y,t)=e iat w(.x,y) 

where w is the frequency of pulsation or oscillation. The 
kernel function K(xo,yo) physically represents the contribu- 
tion to the downwash at a point (x,y) in S due to the presence 
of a pulsating pressure doublet of unit strength located at 
some other point (£,17) in S. It is a function not only of x, 
y, £, and r, but also of Mach number and frequency. The 
function £(£,17) in equation (1) is the unknown lift distribu- 
tion or local doublet strength. (Although it is usually 
convenient to factor out the density term 1/p as indicated 
in eq. (1), this was not done in ref. 1.) 

Equation (1) pertains formally to either subsonic or super- 
sonic flow; however, separate treatments of the two cases are 
required because of wide differences associated with flow 
characteristics. So far as the integral equation is concerned, 
the differences in the two cases he mainly in the kernel func- 
tions. These differences are associated with the differences 
in character of doublets for the two cases. Although the 
main purpose of this analysis is to derive and treat the 
kernel function for the supersonic case, a necessary first step 
is to formulate a doublet suitable for such a treatment. In 
the following section, a desired form, which was arrived at 
by a convenient use of a cutoff, or Heaviside unit function, 
is presented in equation ( 5 ). 

PULSATING DOUBLET MOVING AT SUPERSONIC SPEED 

The governing differential equation for linearized unsteady 
flow at either subsonic or supersonic speeds, which the 
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doublet potentials must satisfy, is the well-known wave 
equation referred to a moving coordinate system: 


bV , dV , bV_I. 

bx 1 by 3 bs 3 c 3 



Under the assumption that disturbances vary harmonically 
with respect to time, this equation becomes 


ay dy sy 1 

bx?~*~by 3 ~^ bz 3 c 3 



where ^ is a complex amplitude function defined by 
Mx,V, z,t)=e ia, $(x,y, z) 


( 2 ) 

(3) 


It may appear that, since the same differential equation (eq. 
(2)) is involved, a logical way of obtaining the potential for 
a pulsating doublet moving at supersonic speed is by simple 
analogy or continuation from the potential for the doublet 
moving at subsonic speed. This procedure is applicable only 
in a broad sense because, as discussed in reference 15 with 
regard to sources in supersonic flow, the potential of a doublet 
moving at supersonic speed consists of the sum of two effects 
corresponding to a retarded-type potential and an advanced- 
type potential which relate to the two wave fronts en- 
countered by a point at any time; whereas for subsonic speed 
only the retarded type of potential is admissible. (The 
advanced-type potential for subsonic or sonic speed does not 
satisfy the Sommerfeld radiation condition, which requires 
that disturbances bo propagated away from their point of 
origin.) In the second place, the potential that may be 
obtained by analogy with the potential for subsonic speed 
must, as subsequently discussed, be rather severely restricted 
before it mathematically describes the physical realities of a 
disturbance moving at supersonic speed. In the following 
development, a desired form of the doublet potential is ar- 
rived at by consideration of these restrictions applied to 
both a retarded and an advanced type of potential that may 
bo obtained by analogy with results for subsonic flow. 

By analogy with results for subsonic speed (for example, 
eq. (A9) of ref. 1) or, more directly, from the discussion of 
source potentials in supersonic flow (ref. 15), the sum of 
the retarded and advanced types of potentials required to 
form the doublet potential for supersonic speeds may be 
written with the doublet situated at the origin as 



-cffi ) 

■ , 

-P 2 * 2 V^-W-iSV 


real values of the radical term are to be con- 

sidered and (b) the values of the expression and its deriva- 
tives are to be considered zero when x is negative. These 
restrictions follow from the physical consideration that 
small disturbances propagate at sonic speed and in a super- 
sonic stream do not progress forward of their point of origin. 

A convenient way of writing the expression for ip D with 
these restrictions accounted for, as previously mentioned, is 
to employ a cutoff or unit function as a factor. Thus, if i f/ 0 
represents the restricted value of i p D , the amplitude function 
of fa may be written as 

fa=2 e- ml *U(x-t sVv^+i 2 ) cos ^ 

where 

U(x—^y 3 +^)= 1 (z>/3- v Y+2 2 )"'| 

, , y (6> 

U(x— /3Vy s +2 2 )=o (z^flVyH-z 2 ) J 

and only positive values of the radical -Jif+z 3 are con- 
sidered. A method whereby this form of potential can be 

determined in a more direct manner is discussed in detail in 
reference 16. The discussion in this reference is in connec- 
tion with the Green functions associated with the dispersion 
of sound waves in an n-dimensional medium in which a 
pulsating source exists. When appropriate changes are 
made in notation, the expression for fa given in equation (5) 
agrees essentially with results for the dispersion of waves in 
a three-dimensional space given in equation (55), chapter 
XVI of reference 16. 

With regard to the unit function TJ{x), in many applica- 
tions where this function is employed it need not be defined 
as having any particular value when its argument is zero. 
In other applications, especially where the unit function is 
involved in a Fourier analysis, it must be defined as having a 
value of % when its argument is zero. In the present case, 
it is conveniently defined, as may be noted in equation (6), 
as having zero value when its argument is zero. 

Derivatives of the unit function give rise to an impulse 
function called the Dirac delta function. For example, 


bx — ^(*)' — 0 (x^O) 

A- U(x)=8(x)= co (x=0) 

A useful integral property of this delta function is 



dx=f( 0) 


o b cos (Af^Var i -flY-/3 2 g 2 ) 


The next step in the analysis is to make use of the doublet 
potential (eq. (5)) to derive the kernel function for supersonic 
speed. 


where M=V/c, /9=-/M 2 — 1, and <o=o>/ F/3 ! . The restric- 
tions that must be placed on this expression are: (a) only 


DERIVATION AND REDUCTION OF KERNEL FUNCTION 

In this section the kernel function is derived and presented. 
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The function is given in terms of an improper integral by 
equation (13) and in a reduced form with no improper 
integrals by equation (15). As it is frequently desirable to 
present results in terms of nondimensional length variables, 
the results given in equation (15) are presented in this manner 
in equation (16). 

In order to derive the kernel function, the function $ of 
equation (3) is considered as the complex amplitude of the 
acceleration potential. As such, \p is directly proportional to 
a perturbation pressure field p=e iat p, through the simple 
relation 

P=— P? (7) 

and to a velocity potential 
through the equation 

y (8) 


By differentiation of equation (8) with respect to z andintegra- 
tion of the result with respect to x, the vertical velocity asso- 
ciated with the acceleration potential \f/ is obtained. Thus, 
when ^ is considered as the potential of a pressure doublet, 
equation (8) affords a straightforward means for obtaining 
an equation for K(xo,y»), namely: 

K(Xii,y 0 )=~ 5(zo,yo,z),_o (9) 


Details of the procedure are as follows: 

The result of the differentiation of equation (8) with respect 
to z may be written as 


y dx dz ^ dz dz 


( 10 ) 


When this equation is considered as an ordinary differential 
equation with dependent variable and independent vari- 
able x, a complete solution is 

{taX foiX 

I = V e ~S-„& (X ’ y ’ s);VdK (11) 

where the lower limit of integration is employed in place of 
a constant of integration and is chosen so as to satisfy the 
condition that ? vanish far ahead of the origin X=0. Thus, 
from equation (9) there is obtained 

j _±5» rx 0 ft _ Q 

K(xo,y o)=Hm ye v J . ° ^ ft>(X,y 0 ,z) e v dh (12) 


or, after substitution of the expression for (eq. (5)) into 
equation (12), the results may be written as 


2 

K(xo,y oHr^lim 
V z ~* 0 




os_(MSVx^^5j) 

yx*-0V* 


dx (13) 


where r= Vyo 2 +z* and, since the integrand is zero for X</9r, 
the lower limit of integration has been changed from — =° to 
P r. 


It is apparent upon examination of equation (13) that, if 
the indicated differentiation under the integral sign is carried 
out, the integrand has singular and perhaps troublesome 
terms. The indicated differentiation with respect to z, 
however, can be replaced by equivalent operations and 
followed by integrations by parts that lead to a reduced 
form of the kernel function containing no improper integrals. 
These steps follow. 

Beduced form of kernel function. — As may be directly 
verified, the indicated differentiation with respect to s in 
equation (13) is, in the limit z-> 0, identical with 


-mx Hm 
2-+0 


U(X — /Sr) cos (A/Zjyx 3 — (HY 1 ) 

yx 2 -j 3V 2 


r 


—l, A 

yjbx V. L yx 2 -dV 

A sin (M5yxw^]}_£^ [§ UtK-PM) 


sin (AfiS-y/x 2 — /S^o 2 ) 


i 

M 


5(X— P\y 0 \) sin (M?yx 2 — 0V)— 


-4AL 5(X — P\yo\) cos (M3^=^V)‘ 
■vX* — P^yo 2 


(14) 


Since the coefficients of 5(X— (% 0 |) in equation (14) vanish 
at X=j% 0 |, it follows from the integral properties of 5(X— /3|j/ 0 |) 
that, when equation (14) is substituted into equation (13), 
the integrals involving the delta function vanish and equation 
(13) becomes: 

J_ e -fc*oU(x 0 —l3\y 0 l) sin (Mu^jxJ— &%*)+ 

~ J^l e-s^Ub-pbo]) sin {M3jxf-p ty 0 *)dxj (15) 

Equation (15) provides an expression for the kernel func- 
tion that involves no improper integrals. Except for the 
integral, the terms of the expression can be quite easily 
evaluated with the aid of trigonometric tables except at 
y 0 = 0, where the function is singular, and at Xo=p\y 0 \, where 
the function is indeterminate. The integral is well behaved 
and can be accurately evaluated by numerical or approximate 
procedures. The singularities and indeterminate values are 
isolated and discussed in a later section, but it is desirable 
first to express the function K(j 0 ,y 0 ) in terms of nondimen- 
sional length variables. As a check on the correctness of 
equation (15), the expression for K(xo,yo ) is reduced to the 
limiting value for M= 1 and compared in appendix A with 
the corresponding limiting value for the subsonic cose. 

The kernel faiction in terms of nondimensional length 
variables. — Although the preceding results contain dimen- 
sional length vadables, it is usually desirable to have such 
results in terms of nondimensional length variables. By 
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employing the variables Xc and y 0 in a new sense to mean that they have been referred to some chosen length l and by 
introducing the reduced-frequency parameter k=fa/V, the length variables may be made nondimensional. (In flutter 
theory the reference length normally is selected as a semichord b.) Tho variables are used in this sense throughout 
the rest of this report. The kernel function (eq. (15)) can be written in terms of these nondimensional variables as 


Itx 0 

K(x 0 ,y 0 )=—y^ e~ ,kz ° [ X ° g C0S (l^ Vzo 2 -£V)-Hg- c ^ U(xo—p\y 0 \) sin V^o*— £W)+ 

ftx 

M £ 0| 6 W) 8 “ 

An equivalent expression for K(x 0 ,y 0 ) which will be useful in subsequent considerations is 

2 


ax o 


(16a) 


K(xo,yo)= 


W 




^tnX TO>r "' % ' l) -(fV^v)]+ 


ax 

U(x ~^ 8in Vx^v) d\ 


(16b) 


Another equivalent expression for K (x 0 ,y 0 ) that is more concise and, for many purposes, more attractive than equation (16a) is 

a _ 


K(x 0 ,y 0 )=- 


Vl'yi 


e -ar ° U(xq— /3|y 0 |) 


... | (io+ Af W-^ a »o s ) 


roe 


cos 




fo-.wVV-flW) 


Vl+r 2 


(16c) 


The reduction of equation (16a) to (16c) is given in appen- 
dix A. One noteworthy feature of equation (16c) is that the 
integral which remains to be evaluated has the same form as 
the integral occurring in the expression for the kernel func- 
tion for subsonic speeds as presented in appendix D of 
reference 1. 

ISOLATION AND DISCUSSION OF SINGULARITIES OF KERNEL FUNCTION 

As previously mentioned and as may be noted in equations 
(15) and (16a), the kernel function becomes singular at 
y 0 =0 and is of an indeterminate nature when Xo=p\y 0 \. It 
is therefore desirable to make special treatment of the 
function in the neighborhood of these values of Xo and y 0 in 
order to be able to express the function in a form which is 
more amenable for calculations. The indeterminate condi- 
tion arises from the first term of equation (16a) because of 
the manner in which the unit function has been defined for 
this analysis. (The denominator of the first term vanishes 
at Xo=/3y 0 . The presence of the unit function in the numera- 
tor, however, renders this singularity indeterminate.) 

In the next few sections the forms of the singularities are 
extracted (see eq. (24)) and tho aforementioned indeterminate 
forms of the kernel function are explicitly determined (see 
eq. (20)). A form of the kernel function more suitable for 
calculation purposes, since the troublesome points are 
isolated, is presented in equation (26). A manner of inte- 
grating the kernel function across its singularities is given 
in equation (27). The singularities of the supersonic and 
subsonic cases are then compared. 

Indeterminate form. — Consideration is first given to the 
indeterminate form, and it is convenient for this purpose to 
consider the value of K{x <i ,y 0 ) at points on the positive 
branch of a hyperbola. (See sketch (b)). 



The equation of the hyperbola may be written as 


ro= « cosh 0 
P\y a \ — £ sinh 0 


} 


(17) 


In these equations £=0 corresponds to x<,=p y 0t since elimina- 
tion of 0 gives x < ) S —p' l y 0 2 =e i . 

After substitution of these expressions for Xq and P\y 0 \ 
into equation (16a), the results may be written as 


K(e,d)= 


2/re 


[2-— it4 COSll 


VI V ginh 2 0 

ike CO?h 5 


f P 

!< 




p 


U(ee~*)e 
Mke 


US mah 6 in 

, , Mke 
cosh 6 cos ■ 


P 


P 2 


sin - 


p 2 


p.cosM _J£ /TUJ. 

A e PU(\~ £ sinh 0) sin( 

J « elnh $ \ P 

■£ 2 sinh 2 0^dxj 


(18) 


To obtain a linpting value of this equation for small values 
of £, the trigonometric and exponential terms can be replaced 
by terms up to the second power of £ in a series expansion. 


43G870 — 07 11 
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If the result of performing this expansion, is denoted by 
K', the equation 


*'M)= 


20 > e~ 1 * t co ’ h e rTt / cosh 9 ik 
~ VP U <« *> IWfiTe' 

k 2 cosh 9 P , cosh 0+l\ 

2p 2fP i0g sinh 9 ) 


(19) 


The first of these integrals may be evaluated from the table 
of Laplace transforms of reference 17 and has the following 
value: 


« * U{\-p\y a \) sin VX 2 -0v) dh=M\y 0 \Kmy a \) 


is obtained, which, in terms of the original coordinates 
ar 0 and y 0 , is 


K'(xo,yo)= 


e -Ux o 

"W 


U(x$- 


Xq 

(2 2ik P\ 

Lv*o 2 -/SV 

W So jSV 


Plo< 


. rp+V^— <3Vi 
jSlj/ol 




(20) 


Although these equations were obtained in order to reveal 
the form of the indeterminate value of the kernel function, 
they are found to contain singularities at y 0 = 0. Prior to 
any further discussion of this result, it is desirable to consider 
the limiting form of K(xn,y 0 ) as y 0 approaches zero, to 
determine all the singularities at y 0 =0. 

Singularities at y 0 =0. — For the purpose of obtaining a 
limiting value of the kernel function for vanishingly small 
values of y 0 , the integral appearing in equation (16a) may 
be written as the sum of two integrals, namely: 


where K x is the modified Bessel function of the second kind, 
of first order. In the second integral, the integrand may 
be replaced by terms up to the second power of y 0 in a series 
expansion. Thus, 


J x e 01 U(\— i3|y 0 |) sin -Jp—py/) d\ 

P" oi i \ / • MkyP Mk\\ ^ 

"* J « 13 Z7(X-j3|yo|) C0S ~fF~) dX 

=U(xo—p\y 0 \) Ji e P (m cos ^ 


Mkxo , • • Mkx 0 \ 
H sm —p- J 



(23) 


XIi e * sin ( 7 ? 

fltx 

= r e~* U(X-p\yo[) sin Vx 2 -0V) d\~ 

Jp\u ol VP / 

e U(k—p\yo\) sin P*yo^ d\ ( 21 ) 


where Oi and Si denote the cosine integral function and sine 
integral function, respectively, which are defined as follows: 



Substituting equations (22) and (23) into equation (16a) 
gives, as a limiting value of K (Xf,,y 0 ), 


lim II(a*,yo)=lum 

po-*0 J?o-»0 


=lim 

tto-*0 


[ Ci (m i) +Ci 

^!7( i ,-*.|){|,+F('r-l) + PlogM- 


(jfe) +< Si (jm )-* Si (afer)]} 
-![ C <3m)+ 0i (^T) +i Kwt)- 4 Si (sfe)]} 


where the following series expression for K x (z) (see ref. 18) is employed: 

z 3 , z 5 , 


^l(z)=(T+l0g|) (-, 


2 ' 16 1 384 


'M-(: 


z . 5z® . 

4 + 64 + 1152" 


) 


(26) 


where 7 is Euler’s constant (7=0.5772157). Examination of equation (24) shows that the only singular terms are the 

same as those which appear in equation (20), namely — and — Pe -lte °log|yo|- Thus, for the purpose of isolating the 

yo 

singularities of the kernel function, only K'(zo,y 0 ), as defined in equation (20), need be considered. Nevertheless, the results 
given in equation (24) may be useful in some applications since they provide a ready means for evaluating the nonsingular 
part of lim K (xo,yo)- 

yo-*0 
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Form of kernel function suitable for calculations. — As in the subsonic case, with knowledge of the critical values of 
the kernel function, an expression can be written in which the kernel function is separated into two parts, one of which 
contains no singularities or indeterminate values and the other of which contains all the singularities and critical values 
of the kernel function. This expression is 

K(xo,y 0 ) = [Kixo.yd—K'ixoMoJl+K'ixoMo) (26) 

where K(xo,y 0 ) is defined in equations (15) or (16), and K' (xa,y 0 ) is defined in equation (20). The term K (xo, y Q ) —K' (xo,yo) 
in equation (26) has no singular or indeterminate values. The term K'(xo,y 0 ) is singular at y 0 =0 and indeterminate when 
%Q=pyo. 

Integration of singularities of kernel function. — Since integration of the kernel function is often necessary, a few remarks 
on how to circumvent its inherent singularities are in order. Each term of K' (to,2/o) in equation (20) possesses a simple 
indefinite integral with respect to the variable y=y—yo- In performing integrations with respect to t] that involve a 
passage across the line y =y, a principal value is to be taken. Eor example, 


U (to— P\yo\)x<> drt l 

U (x$ p\y 1 1) sjx^—p 2 (y— l) 2 U (to— 0|y+l |) V^o 2 — d 2 ( 3 /+I ) 2 

J-i (y— rif-jxj—p 2 (y— j/) 2 *o 

y — i 2 /+i 


(27) 


where the symbol indicates that the singular integral is to 

be considered simply as a function of its limits. A justifica- 
tion for this consideration is that it leads to results that 
could, with considerable labor, be rigorously established by 
maintaining the variable z in the analysis until all operations 
aro performed. 

Comparison with singularities of subsonic case. — It may 
be of interest to compare the above results with correspond- 
ing results for the subsonic case, that is, for To>0 and 
j/o »0. Results for the subsonic case may be obtained 
from equation (31) of reference 1 as follows: 


ft 

-ay-MVa^-Hl— M*)y 0 a 


ik 


Vzo 3 +(l-M 3 )yo 2 2 (1 —M 1 ) ^+(x-M t )y 0 i 

ft* , . _ MW+ {i-my^-x,) \ 

2 l0g 2 0=34) J 

r 2 ik J<? P , k(l+M)y o n 
~VP L yo 2(1 +M) 2 l0g 4ro J (28) 

The singular terms of this expression for Xo|>0 and M<( 1 are 


of the reduced-frequency parameter (see eq. (29)) and an 
expansion in terms of Bessel functions. The latter expan- 
sion is used in a later section to obtain the kernel function 
for two-dimensional flow from that for three-dimensional 
flow. 

Power-series expansion with respect to k . — As in the case 
of subsonic flow, the kernel function can be expanded into 
a power series with respect to k that, in the present case, 
is useful for small values of k/p 2 , a combination of reduced 
frequency and Mach number that is inherent in such an 
expansion of the supersonic kernel. The terms of the expan- 
sion may be simply obtained by expanding the terms of 
equation (16a) that are functions of k and collecting the 
results. The first few terms are 


* (*o,yo)=-|2 U(xo—p\y 0 \) 


Vl*y 0 ~ 


ik fPyf 


1 fk 


fiVo 2 0*V*o S — W 

5 |)‘ 


-i - a i } - \+ 


COsl1 

1 f AY (3M i -l)dVoV+(2-3M 2 )<3V , n Y 
6VJSV Va--o a -/3V + W. 


(29) 


Comparison of this result with equation (20) shows that the 
singularities for subsonic and supersonic flow are of identical 
form. 

SOME INFINITE-SERIES EXPANSIONS PERTINENT TO THE KERNEL 
FUNCTION 

The kernel function can be expressed as a series by various 
expansion procedures. Some particular expansions, which 
should be useful in applications, are discussed in succeeding 
paragraphs. These are the power-series expansion in terms 


Although this power-series expansion converges to the 
appropriate value of K(xo,y 0 ) for all finite values of k/p 2 , a 
great number of terms are required unless k/p 2 is small. 
These first few terms of the expansion can be considered to 
represent the kernel function for values of k in the range of 
magnitudes generally encountered in dynamic-stability 
studies and, therefore, they are pertinent for obtaining 
time-dependent stability derivatives. A noteworthy feature 
of the expansion is that each term can be integrated, in the 
sense that it contains a simple indefinite integral, with re- 
spect to the variable i]=y—yo- When such integrations 
involve a passage across the line y=y, a principal value is 
to be taken in the sense described after equation (28) . 
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Expansions in terms of Bessel functions. — The trigono- 
metric terms appearing in the expression for K(xo,y 0 ) in 
equations (15), (16a), and (16b) can be expanded into infi- 
nite series involving Bessel functions of the first kind. Such 
expansions have good convergence properties, even for large 
values of the parameter left S 2 , and each term possesses a 
simple indefinite integral with respect to q. Such Bessel 
function series are therefore useful for deriving an expansion 
for the indefinite integral of K(xo,y 0 ) with respect to q. 
The indefinite integral of K(xo,y 0 ) leads to the downwash 
associated with pulsating vortex lines (“horseshoe” vortices) 
and, as previously indicated, to the kernel function for two- 
dimensional flow. It might be useful to point out that the 
expansion of the cosine term into a series involving Bessel 
functions is also useful for studying distributions of pulsating 
sources. 

For the purpose of expanding the trigonometric terms 
under discussion, consider the expressions 

Cf(X-„)S21 ME? (30) 

and 

U(\—a) sin (&yx 2 -a 2 )=— Jj? 1 (X— a) 

whore a, b, and X are positive. 

By making use of a known Fourier transform relation, 
expression (30) can be equated to an infinite integral involv- 
ing a Bessel function of the first kind (see, for example, 
p. 33 of ref. 17) : 

f“ Jo (xVAfP) cos ardr^ C( ^ySE^ 

Jo yx 2 — ar 

J Jo (xV?T^) COS ardr= 0 


(X>o) 

(x<a) 


Y (32) 


By use of the addition formula for Bessel functions (see. 
for example, p. 358 of ref. 18), the Bessel function appearing 
in this equation can be written as an infinite sum of products 
of Bessel functions as follows: 

JqCxVt 2 +i 2 )=t7o(rX)Jo(6X)-4-2y^, ( — l) n JiJj\)J2„(b'K) (33) 

n-1 

Thus, 

P Jo (XV^+F) cos aTdr=U(\-a)^^^^^- 
J o y X 2 — a 2 

~J 0 [Jo(tX)Jo(6X)+ 

2 2X — 1)V 2b (tX)J Sb ( 6X)1 cos ardr 
n-l J 

(34) 

In view of the relation (see ref. 17, p. 37) 

U(\—a) cos fin sin -1 
JsJrX) COS ardr== 

Jo yx*— a . 2 

the indicated integration on the right-hand side of equation 
(34) can be carried out term by term so that 


U(\—a) cos (fr-y/X 2 — a 2 ) 

Vx 2 — a 2 

I 2 [' 7 °( 6x )+ 2 S (— 1 )*Jj»(JX) cos (2n sin" 1 (36) 

Substituting the expression on the right-hand side of equation 
(36) into equation (31) gives 


U(\—d) sin (6VX 2 — Ji(b\)— jb (— 1)" Jj«-i(6X)— Jj* +i ( 6X) cos^2n sin -1 


2XZ7 (X— 
’ V>r 


cos (sin *0S(— 1)* 'JiM-iibX) cosj^(2n— 1) sin'^JJ. (37) 


But, since 


U(\—d) ^j yi _ a 2 cos (bin -1 ^U(\~a) (a^ 0) 

the expression for U(\—a) sin (jy\ 2 — <z 2 ) may be written as 


U(\~d) sin ( b VX 511 ^) =2Z7(X — a) Zj (-l)*" 1 Jt»-i(.b\) cos |J2»- 1) sin' 1 £] (38) 

By direct comparison, equations (36) and (38) can be used to write expanded forms of the trigonometric terms appearing 
in equations (15), (16a), and (16b). Expansions thus obtained will now be used to derive the kernel function for two- 
dimensional supersonic flow. 
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In contrast to three- dim ensional flow, a physical interpretation of the kernel function for two-dimensional flow is that 
it represents the downwash at a given field point due to a pulsating bound vortex line of infinite length. This kernel 
function may be obtained by integrating the kernel function for three-dimensional flow from— <=to or, in view of the role 
of the unit function, from one Mach line to the other, with respect to the variable v=y—yo ■ Pulsating “horseshoe” 
vortices may be obtained by integrating K {x a ,y Q ) over an arbitrarily finite range with respect to i?. 

DERIVATION OF KERNEL FUNCTION FOR TWO-DIMENSIONAL FLOW 

In this section the kernel function for three-dimensional flow is reduced to the function for two-dimensional flow, and the 
final results of the reduction are given in equation (49). For the purpose of derivation, K(xo,yo ) will be considered as given 
in equation (16b). The two-dimensional kernel function can then be expressed as 


P 1 


J“ K^yo^-yi | « [tf(*o-/SN) sin (j~ 

b $ L U{ ^ v « sta (f v?=w ) 



=K(x 0 ) 

(39) 

or 

2 / k \ 




(40) 

where 




\utz«-p\ y -r>\) sin V^-/?Ty-n) 2 )] dv 

(41a) 

or 

(m l+n) £ *“ c?^) * 

(41b) 

and 

h=P ^ Qj* « ** U(K-\ V |) sin V^=V) dx] dv 

(42) 


In equations (41) for /j and equation (42) for I 3) infinite 
limits with respect to the integration of v are not necessary. 
In view of the role of the unit function in the integrands, 
limits of ±Xo include all values of v for which the integrands 

are different from zero. The symbol ^ indicates that the 

integrals are to be considered simply as functions of their 
limits or that the singularity at 17=0 is to be ignored. 

First consider equations (41) for J x and then perform an 
integration by parts. The expression for /j may then be 
written as 


By direct comparison with equation (36), equation (43) may 
be written in expanded form as 



2S(-1)”4 

n— 1 


(Mk 

VF* 


cos (2n sin -1 dv (44) 


In this equation, the terms involving J 2n do not contribute 
to the integral because 


l*b 

-*0 


m ss+b) | 

— ! 2 K^ + f)X"^? ' (43) 


f T ° cos (271 sin -1 — y dr] 

Jo — rf \ *0/ 

j* 0 j^jj 8(xo— ]u|) sin ^2n.sin -1 di/JsO 
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f r ° d v =Ufa-\v\) sin -1 ^ *+ 

Jo 11 *0 0 

tlie expression for I x can be 'written as 

=—ttJ 3 [« (atf +|f Ufa) Jo -Mb Ufa) J x *,)] 

(45) 

Now consider equation (42) for I a , namely 


l — i£X 

L=P fi 7 |_Jm ^ 5111 (f v^=?) <&] ^ 

The double integral in this expression can be considered as a 
surface integral over a triangular region of the X 77 -plane cut 
out by the lines v = \ v—~ \ and \=Xo, as shown in 
sketch (c). 



Sketch (c) 

By a change in the order of integration, which is admissible 
since the singularity at 77 = 0 is to be ignored, the expression 
for Ja may be written as 

I 3 =/3 X" 0 ^ ^11^7^ 8hl d v (46) 

The inner integral in this equation is identical in form to the 
integral in equation (41). Hence, by observation of and 
comparison with the results obtained for _Zx in equations (43), 
(44), and (45), it is found that 



The expression for I a can therefore be written 

la= ~ ^ Ufa) Jo x) dk (48) 


Substituting this result and the results given in equation 
(45) for Ji into equation (40) gives a desired form of the 
kernel function for two-dimensional supersonic flow : 



Examination of equation (49) shows that the only singularity 
involved in the kernel function for two-dimensional super- 
sonic flow is the 5-function. At zero frequency, all the 
terms of Kfa) except the 5-function vanish. The kernel 
function required to treat two-dimensional wings at steady 
angle-of-attack conditions is therefore proportional to this 
5-function, and, as shown in the following section, leads in 
a very simple manner to the well-known Ackeret results. 

The integral that remains to be evaluated in equation 
( 49 ) is well behaved and similar to integrals, treated by 
Schwarz (ref. 19) and others, that arise in the velooity- 
potential approach for treating two-dimensional wings. 

APPLICATION OF KERNEL FUNCTION TO LIFT DISTRIBUTIONS FOR TWO- 
DIMENSIONAL WINGS 

The results obtained in the previous section for the two- 
dimensional kernel function are now employed to obtain the 
lift distribution on oscillating and steady two-dimensional 
wings moving at supersonic speed. (See eqs. (56) and (61).) 
Since the lift distributions so obtained agree with the Ackeret 
results for a steady wing and also -with known results for the 
oscillating wing (ref. 15), they serve as a check on the correct- 
ness of the expressions for both (ho two-dimensional and 
three-dimensional kernel functions. 

The integral equation that must be solved to obtain the 
lift distribution for two-dimensional wings in supersonic 
flow is particularly simple since it involves a single integral 
of the convolution type: 

w(xfa-^- p £ mt)Kfa)dt=±- p j Z o L(£)K(x-m ( 60 ) 

Integral equations of this type can be readily solved by 
Laplace transform procedures since the Laplace transform 
of a convolution integral is the product of the transforms 
of the functions that compose the integrand. In the present 
case, if 8 represents the Laplace transform operator defined 

*>y 

£L/(z)]=J o “ e-*f(x)dx=M (51) 

the transform of equation (50) may be written as 

w{8)=~U8)K{ 8 ) (62) 

Solving this equation for L(s) gives the Laplace transform 
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of the lif ^distribution : 


^- 4 -F 1 


Inversion of the transform on the right-hand side of this 
equation gives the lift distribution. 

For the^caso of a steady two-dimensional wing, 


w(x)=Va 


K(x)=^f 5 (x) 


w(s)=— 

s 


m=- 


L(s)= 


2 P V' 2 a 


The inveme transform of equation (54) gives for the lift 
distribution: 

Ux)=^^U(x) (55) 

From this result, the total lift per unit of span is 


\sliui u 

ij L{x)dx— 


2pV 2 aX Chord 
P 


This result agrees with the well-known Ackeret result. 

Now consider the unsteady case for oscillatory translation, 

w(x)=ialh=iVkh. (57) 

where h is the amplitude of displacement referred to Z, and 
the Laplace transform of w(x) is 


, iVkh 

W(8)= 

S 


The Laplace transforms of the different terms of K(x) 
(eq. (49)) can be simply derived or they may be obtained 
from Laplace transform tables (for example, ref. 16). After 
combining the transforms of the different terms, the results 
can be written as 



Substituting equations (58) and (59) into equation (53) gives 


for the transform of the lift distribution 

Us) 


The inverse of this transform gives for the lift distribution 

<6i) 

This result can easily be shown to check with the results of 
reference 15. Moreover, if ikh is set equal to a, and then k 
is allowed to approach zero, equation (61) reduces to the 
result for the steady case. 

CONCLUDING REMARKS 

The main purpose of this report was to derive and present 
in a form that could be numerically evaluated the kernel 
function of the integral equation relating downwash and lift 
distributions for oscillating wings in supersonic flow. This 
purpose has been achieved for three-dimensional flow, and 
the results have been converted to a form more suitable for 
calculation by isolating the singular or critical points. The 
kernel function for two-dimensional supersonic flow has been 
presented and the results show that the only singularity is a 
Dirac delta function, which appeared in such a manner that 
further reduction with regard to singularities is not required. 

The results presented in this report for supersonic flow 
together with those previously obtained for subsonic flow 
provide a kernel function that is capable of being evaluated 
at any Mach number. As experience develops it is expected 
that use can be made of the kernel function to develop 
approximate procedures, that will be more or less uniform 
throughout the Mach number range, for calculating aero- 
dynamic forces on oscillating (or steady) wings of arbitrary 
plan form and with arbitrary downwash conditions. The 
labor involved in such approximate or numerical procedures 
will indeed be prodigious and will require the use of modem 
high-speed computing equipment. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 15, 1955. 


2ipV i kh 

P 


// , iMW 



158 


REPORT 1257 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX A 

DERIVATION OF EQUATION (16c) AND REDUCTION OF THE KERNEL FUNCTION TO THE SONIC CASE 


The purpose of this appendix is first to reduce equation 
(16a) to equation (16c) and then to reduce equation (16c) 
to the sonic case. The reduction to the sonic case, by com- 
parison with results obtained from consideration of subsonic 
speeds for the sonic case in reference 1, provides a partial 
check on the correctness of results obtained for the supersonic 
case. 

DERIVATION OF EQUATION (16c) 

In order to reduce equation (16a) to equation (16c) con- 
sider the integral 

^ Xj ° i e ° r ' sin Vx 2 — 7 x 

(Al) 


-3 

where 


and 


I t = ( X ° U(\-p\y 0 \)e dx (A2) 

I a = P U(\~p\y 0 \)e~^ x+M ^ r ^ ) dk (A3) 

In these expressions X takes on only positive values, that is, 
P |y 0 | ^X^xo. Hence, consider for each integral the single- 
valued substitution 


X= lyoK+AfVl+T 3 — t) 

which, for leads only to positive values of X. 

Solving this expression for r gives 


(A4) 


X±My/X 2 — 0 2 y o 2 

T ~ W 

Thus, if the substitution 

(X— j3Vo s )=/3% 0 ]r 
is made in Ii and the substitution 

(X+My/X 3 — /3 s 2/o 2 )=/3 ! ji/o|r 


(A6) 


(A6) 


(A7) 


A(xo,y 0 )= 


vl s yo' 


e % 0 |) 


is made in h, there is obtained 
,^j (xo-MW-W) 






and 
Ii 


(A8) 


=wj‘ 


■^1 




u HXw; i- 1 ) 8 -" 1 -"*' 


= W J. ir+ 


VVl+T 7 


T dr 

(AO) 


tC 

Combining these results gives for (7i— Ii) 
2^(Ii-/ 2 )= 


MyolUixo-PM) p’W (ro+WV_W) / Mr 

2iM J , . __ , VATT 3 / 

^j-, (xo-mV^w) 7 

(A10) 


or since 


/’flliTl O r o+ M vV-/ s W) 

r w e - tt i»okd T= 

(xo-AfV^ 3 -^ 1 ) 


2§ ? (/ l -/ 2 )=-yC7(xo-/3h/oD 


2e . IMk 


sin 0^ V^o 2 — /3V^) 
(All) 


it 


ie *’° . (Mk 


^V*o s -/3W) 




+ 


(A12) 


Substituting this result into equation (16a) gives equation 
(16c) or 






V* o 2 — /3*2 /o : 




cos 




(ij-Afy/V-PW) 


Vi+ 7 


, g— tt k«l T d r 


(16c) 
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REDUCTION OF EQUATION (16c) TO THE SONIC CASE 


In order to reduce this result to the sonic case it is first necessary to discard terms arising from the advanced-type 
potentials employed in deriving the doublet potentials for supersonic flow. This may be accomplished by replacing 

/Mk i—5 A . l 

cos {-jp ^XQ—^yaj with ^ & • The limiting form of K(xo,yo ) as M— >1 can then be written as 


^o,yo)j/~i=lim— -T5- 5 e~ itx o U(xo—p\y 0 \) 
Vv y o 


Xq ~jji (io -M-yJxj-pij) 

e 




r 2 M 


iTTi (xo+l-fv'^W) 




-JT+t* 


e-'Voh dr 


(A13) 


When the limit M= 1 is approached from the supersonic side, the term M is conveniently replaced by 

M=l+|t 

where < is infinitesimally small, so that 

/S 2 =(M-l)(M+l)=i e (2+i 

With this approximation equation (A13) can be written as 


(A 14) 


K(x o,yok-i=lim-^- 2 U(xo-^\y 0 \)-{ 


xo 


■\]xo — ey 0 i 


P itx o 

=lim — 5 U(xo) 

w, vPy 0 2 




t[-o- (!+-;) vv^ 5 ] r ih + ( l+ 0vv-^] 

+*%o|J 

^ ~ ( 1+ i) V^o* - ] 


Vl + r 2 


dr 


VT+^ 


e <i i*'ol T dr 


(A15) 


or 






+i% 0 |f u N r dr 

j 1 ( vo*\ yi+r 2 

V 1 ” *0 / 


(A16) 


In order to show that this result is equivalent to that given for the sonic case in reference 1, it is convenient to first 
express the integral term as the difference of two integrals 


/% 03 

*%o|J 


~%I 


H 9 


-jl.citww/, 




e _,i I* , ol r dr- 


-U\yo\( 




Vl + r 2 


e- tt ^oh dr (A17) 


The first integral on the right of equation (A17) can be evaluated by comparison with the following results: 

J '° f r® j # f*to ^ 

■ ; e~ tar dT= ■ , cos ar dr — i , — = sin ar dr 

0 ^it+? jo jo vr+r* 

-|Y r^dr+i r^*) 

^°\Jo -\/l + t 2 Jo -JT+r 2 / 

= da"^2 t A(a-) — Ao(a)] + iK^a) J- 


=2 [A(a) — ii(c)] — 1 —iK^a) 


(A18) 


435875—57 12 
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where Io, Iu Ko, Kr are modified Bessel functions and Lq and Li are modified Struve functions (see ref. 18, p. 172 and p. 
332 for reductions of similar integrals). Thus 

ttfe/o| J 0 " e~^o\r dT ^y 0 \ Kl(Jc \y Q \ )+ ^M [7 1 (% 0 |)-X 1 (% 0 |)]-i% 0 | (Al0) 

In the second integral on the right of equation (A17) make the substitution 



or 

This gives for the second integral 


^=ll/o| G/l+r 2 — t) 



(A20) 


Substituting this result and that given in equation (A19) for the integral in equation (A16) gives 

^ U(x a )j\2et [ I «^]_i-iij 2/o | + ij yo jZ- 1 (^|y4) + ^|^ [imo\)-In{k\yo\)]-ik^^ J ( X_ ^) j (A21) 

A comparison of this result with the result given in equation (47a) of reference 1 shows that the two equations are 
equivalent. 


APPENDIX B 

DERIVATION OF DOWNWASH FUNCTIONS ASSOCIATED WITH “HORSESHOE” VORTICES IN SUPERSONIC FLOW 


The downwash associated with a vortex line can be 
obtained by an integration, between appropriate limits, of 
the kernel function K(x 0 ,yo) with respect to ij=y—y 0 . In 
order to perform such an integration analytically, recourse 
must be had to term-by-term integrations of an expanded 
form of K(xo,y 0 ). In this regard, use can be made of the 
expansions given in equations (36) and (38) of the analysis to 
obtain expanded forms of the downwash functions for vortex 
lines that have very good convergence properties, especially 
for the range of values of the parameter Mkj\ 3 1 that would 
usually be of interest in applications. Expressions so ob- 
tained will be cumbersome and will require high-speed 
computing equipment to make them very useful. 

In regard to “horseshoe” vortices in supersonic flow, there 
are five different significant regions in which a field point 
may be considered to be located (see sketch (d)). 

Region (1) is between the Mach cones emanating from the - 
end points of the bound-vortex line. The downwash at a 
point in this region is not affected by the trailing vortices but 
is created by the bound vortex alone. Therefore, the down- 
wash is the same as would be produced by a bound vortex of 
infinite length and corresponds to the kernel function for 
two-dimensional flow discussed in the analysis. Region (2) 
is between the trailing-vortex lines and is within the Mach 
cone emanating from one end of the bound vortex but outside 
the Mach cone emanating from the other end. The down- 
wash at a point in this region is created by the bound vortex 
and one of the trailing-vortex lines. The other trailing 


( 4 ) 


Trailing vortex 


/ 

/ ( 5 ) 


( 2 ) 


Wind direction 


(I) (3) 


( 2 ) 


Trailing vortex (5) 


( 4 ) 


Sketch (d) 


vortex has no effect on the downwash. Region (3) is be- 
tween the trailing-vortex lines and is within the Mach cones 
emanating from both ends of the bound vortex. Downwash 
in this region is created by the bound vortex and both trail- 
ing-vortex lines. Region (4) is outside the trailing-vortex 
lines and is within the Mach cone emanating from one end of 
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the bound vortex. The downwash is created by the bound vortex and only one of the trailing-vortex lines. Region (5) is 
outside the trailing-vortex line but within the Mach cones emanating from both ends of the bound vortex. The downwash 
is created by the bound vortex and both trailing-vortex lines. 

For any of the five regions discussed in the preceding paragraph, the integral corresponding to the downwash function 
may be formally written, with use of equation(16b), as 

r *** „ jia+jl 

jj 1 K( ^>Vo) dr l=~y^ e ^ ffk/o|) sin VV~ ffW) dv+ 

ux “I 

CM) 

where uso of the substitution j3j/ 0 =l gives 



In equations (B2) and (B3), a principal part — as described after equation (27) in the analysis — is to be taken when the 
integrations are carried across the line f=0. The purpose now is to reduce these expressions to forms amenable to 
numerical evaluation. The first step in this procedure is a reduction of the expression for I 4 (eq. (B3)). The double 
integral in this expression can be considered as a surface integral in the Xf-plane where the order of integration is first with 
respect to X and then with respect to f . The steps in the reduction are first to delineate the area of integration for each of 
the five different cases under consideration, and then to change the order of integration in the surface-integral represent- 
ation of It. 

From the description of the different cases to be considered and by examination of the limits of integration in equations 
(Bl), (B2), and (B3), the area of integration for the case of a field point in each of the aforementioned regions may be 
considered as shown by the hatched areas in the following sketches: 



( 3 ) ( 4 ) ( 5 ) 

Sketch (e) 
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Expressions for /« for the five different regions or cases may then be expressed as simple integrals as follows: 
Case (1): 

r* — 

I 4 =J e ? 


Case (2): 

Case (3): 

Case (4): 

Case (5): 

The expression 


F=£'e 01 * F(\,-^)dK 

rh rt , r*o — 

I 4 =J o e * F(X ,-\)dk+ s Ffo-fOdX+J^ e F(r 2 ,-ft)dX 

I *= 6 0 F(\,£i)d\ 

Tfs — C x — — 

F(X,fi)dX+ ^ FfaJJdX 

F(\,-\)= sm vx’-r* * 


(B4) 

(BS) 

(B0) 

(B7) 

(B8) 

(BO) 


is evaluated in the text in connection with the derivative of the two-dimensional kernel function and is found to reduce 
(see eq. (47)) to 

F(\-\)=^U(\)Jo(~: 


0 


(BlO) 


The F-function for other arguments can be obtained by substituting appropriate limits in an integration by parts of 
F(X,—f i), namely (see the development following eq. (41) in the analysis): 

J-s)? ^ 1 1|} °in Vg=? it 

— ]OTHrl) sinf: J.(^) TOHrl) A (nr), lj! r ,ifj »<x-lfl)»in- [<K- 

^ A. (^) [w-lrl) <*■ (2» la- £) f_ t + J^jf, 5(Mfl) sin ( 2 * sin- £) if] (Bll) 

After the first term on the right-hand side of equation (Bll) has been expanded by comparison with the expansion 
given in equation (38) of the analysis and the limits of integration have been substituted, this expression may bo 
written as 


F{X,~U)=^ U(\) Jo (^)-Z7(X— |M) Jo (~) sin- 1 £+S [ 2( Jin-\ (jpr) cos [(2a-l) sin" 1 £]+ 

Mk (-1)" T /AftX\ . f n . 

sm V 271 sm x)J/ (® 12 ) 

Substitution of f 2 for X in the limits of equation (Bll) gives 


F& ,-f 1 )=-Z7(X-|f i |)^ Jo (~) sin" 1 f +g [ 2( ~ r f _1 (^) cos [(271-1) sin" 1 ^]-!- 

ds- (^pn) sm (2 ti sin" 1 §)] }~-O r 0HiiD-^jr do (jyr) sm" 1 y+ 

- Jin-i(—^r ) cos [(2n— 1) sin- 1 ^]-^ j in sin ( 2 n sin- 1 £)]} 


.Mfc (— l) a 
/3 2 71 

'2(— l)"- 1 


(B13) 
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If — ft is replaced by ft in the limits of equation (B 1 1 ) , then 


F(Ktd=^U(\)J 0 


/Mkh\ , 

V p 2 r 




sin" 1 £+ 
K 




'2(— 1) 


n-l L fi 


Jin - 1 cos £(2n— 1) sin -1 yj+ 

Mk (-l) n T / Mk\\ . /_ . 

^ sm ( 2nsm yj/ (Bi4) 

Substitution of ft for X and ft for — ft in equation (Bll) gives 


Fit*, ft)=-?7(X-|ft|)^ Jo (^) sin- 1 £+ 

§ [ 2( ~ft J ' 2 - 1 (^) 008 [^-i) sin_1 $]+ 

J ’* (*f) sin ( 2n 8in_1 r)]} + 

^-lr 1 |){^^o(^)sin- i ^+ 

§ Fir^ Js - 1 (nr) cos [ (27l ~ 1} sin “ 1 f] + 

M/ ^n 1) ' 1 J ’" (^) 8in (^sin- 1 £)]} (BIS) 


When equation (BlO) and equations (B12) to (Bl5) are 
substituted into equations (B4) to (B8), respectively, they 
give the reduced forms of I t for the five cases under considera- 
tion. 

After the reduction of It, the corresponding reduction of 
I 3 is considered. As may be found by examination of the 
expression for I 3 (eq. (B2)) and the sketches showing the 
areas of integration for the different cases, reductions of 
I 3 corresponding to those of can be obtained from the 
A-functions (eqs. (BlO) to (B15)). Results for the different 
coses may be expressed as follows: 

Case (1): 

/ 3 = 1 ™ (wc W*ll) ' F(X) _X) (Bl6) 

Case (2): 

/3= ££ (mi $k + M) c 817 ) 

Cose (3): 

/3= }i“ (mc &x + m) (Bl8) 

Case (4): 

® 19> 

Cose (5): 

® 2I » 


When the expressions for I 3 (eqs. (Bl6) to (B20)) and I 4 
(eqs. (B4) to (B8)) that axe associated with each particular 
case are substituted into equation (Bl), expressions for the 
downwash at each of the five significant field-point locations 
may be obtained in terms of the A-f unctions (eqs. (BlO) to 
(B15)) as follows: 

Case (1): 

X. K ^y^ di > = —f\_ e * S? # (m ik+m) F(k ~^+ 

— i 

M e itXJ J 0 e ** A(X,— X)dxJ 


(B21) 


Case (2): 


lAfSfcn, 


J: K^,y 0)dv= -fh- » Hm(&|4)^-ft) + 


®y/: 


0 

ikX 


(B22) 


1 e * F(X,~X)dX+ 

J f y^^(X,-ft)dX 

Case (3): 

f iM*kXo 

L 01 £? # (m w&m) F fa-tt+ 

^[j^F^ + 

(IX 

JS“^A(X,-ft)dX+ 


ux 

J*e F(t a ,-fOd\ 


Case (4): 

K(xo,y 0 )dT]= 
Jv 1 


M 

i ' 


iJVUfa,, 


X-.x„ VI 


Mk b\^M 

fjix 


(IX 


(B23) 

y F( j’ fo+ 

(B24) 


Case (5): 

J ”‘ K(xo,y 0 )dv= 

Vl 


|= — e pl lim 

t X-*r, 




Jl 

M 


e~ Ux o 



ikX 

e F(\ ft)dX+ 



(IX 

e ^ A(ft,ft)dX 


} 


(B25) 


The results for case (1) (eq. (B21)) agree with results obtained 
for the two-dimensional kernel function given in equation 
(49) of the analysis. 
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